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Abstract 



The conditions for the existence of negative-energy perturbations (which could 
be nonlinearly unstable and cause anomalous transport) are investigated in the 
framework of linearized collisionless Maxwell-drift kinetic theory for the case of 
equilibria of magnetically confined, circularly cylindrical plasmas and vanish- 
ing initial field perturbations. For wave vectors with a non-vanishing compo- 
nent parallel to the magnetic field, the plane equilibrium conditions (derived by 
Throumoulopoulos and Pfirsch [Phys Rev. E 49, 3290 (1994)]) are shown to re- 
main valid, while the condition for perpendicular perturbations (which are found 
to be the most important modes) is modified. Consequently, besides the toka- 

mak equilibrium regime in which the existence of negative-energy perturbations 

d In T 

is related to the threshold value of 2/3 of the quantity rj u = ^j^jjyS a new 
regime appears, not present in plane equilibria, in which negative-energy pertur- 
bations exist for any value of rj u . For various analytic cold-ion tokamak equilibria 
a substantial fraction of thermal electrons are associated with negative-energy 
perturbations (active particles). In particular, for linearly stable equilibria of a 
paramagnetic plasma with flat electron temperature profile (r) e = 0), the entire 
velocity space is occupied by active electrons. The part of the velocity space occu- 
pied by active particles increases from the center to the plasma edge and is larger 
in a paramagnetic plasma than in a diamagnetic plasma with the same pressure 
profile. It is also shown that, unlike in plane equilibria, negative-energy perturba- 
tions exist in force-free reversed-field pinch equilibria with a substantial fraction 
of active particles. The present results, in particular the fact that a threshold 
value of i] u is not necessary for the existence of negative-energy perturbations, 
enhance even more the relevance of these modes. 



I. INTRODUCTION 



Negative-energy perturbations are potentially dangerous because they may be- 
come nonlinearly unstable and cause anomalous transport - [[HJ. Conditions 
for the existence of perturbations of this kind can be obtained on the basis of 
the expressions for the second variation of the free energy which were derived by 
Pfirsch and Morrison |J for arbitrary perturbations of general equilibria within 
the framework of collisionless Maxwell- Vlasov and Maxwell-drift kinetic theories. 

For homogeneous, magnetized plasmas and vanishing initial field perturba- 
tions they found that negative-energy perturbations exist for any wave vector k 
having a non-vanishing component parallel to the magnetic field (parallel and 
oblique modes) whenever the condition 



>0 (T 



dv\\ 

holds for the equilibrium guiding center distribution function fj® for some par- 
ticle species v and parallel velocity v\\ in the frame of lowest equilibrium energy. 
For inhomogeneous magnetically confined plasmas with equilibria depending on 
just one Cartesian coordinate y, Throumoulopoulos and Pfirsch showed that, 
in addition to parallel and oblique modes, for which condition (1) applies, per- 
pendicular modes have also negative energies if 



dP<°> dfj® 



< 0, (2) 



dy dy 

holds, where is the equilibrium plasma pressure. For tokamaklike equilib- 

d In T 

ria, condition (2) implies a threshold value of 2/3 of the quantity rj v = Q\ n pf i 
where T v is the temperature and N v the density of particle species v. These 
investigations are extended in this paper to the more interesting case of circu- 
larly cylindrical plasmas. The method of investigation consists in evaluating the 
general expression for the second-order perturbation energy obtained by Pfirsch 
and Morrison within the framework of the linearized collisionless Maxwell-drift 
kinetic theory. The most important conclusions are: 

1. Condition (1) for the existence of parallel and oblique modes remains valid. 

2. For tokamak and reversed-field pinch cold- ion equilibria a new regime ap- 
pears, not present in plane equilibria, in which perpendicular negative- 
energy perturbations exist without restriction on the values of r\ v . 
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The equilibrium properties of the circularly cylindrical plasmas under con- 
sideration are discussed in Sec. II. The second-order perturbation energy for 
vanishing initial field perturbations is presented in Sec. III. The relevant lengthy 
derivation is reported in Appendix A. The conditions for the existence of negative- 
energy perturbations are obtained in Sec. IV. The cases of parallel, oblique and 
perpendicular wave propagation are examined separately. The consequences of 
the condition for the existence of perpendicular negative-energy perturbations in 
straight tokamak and reversed-field pinch equilibria are discussed in Sec. V. For 
various analytic cold- ion equilibria with non- negative and negative values of r] e , 
the part of the velocity space occupied by electrons associated with negative- 
energy perturbations is also obtained. Two examples are presented in Appendix 
B. The main results are summarized in Sec. VI. 



II. EQUILIBRIUM 



The collisionless Maxwell-drift kinetic theory applied in the present paper is 
based on Littlejohn's Lagrangian formulation of the guiding center theory |16| in 
the form given by Wimmel |Tj]]. A brief review of this theory is given in the first 
paragraph of Sec. III. More details can be found in Ref. and in Sec. II of Ref. 

For a magnetically confined, circularly cylindrical plasma the equilibrium vec- 
tor potential and magnetic field are given by 

A(°) = 4 0) (r)e e + 4°)(r)e, (3) 

and 

B(°) = 5f (r)e e + 4°)(rK, (4) 

with 

-UrAf) = B?\ {Af)' = - B f\ (5) 
r ar 

Here, r, 8, z are cylindrical coordinates with unit base vectors e r , eg, e z and 
the prime denotes differentiation with respect to r. It is assumed that there is 
no equilibrium electric field. To calculate the guiding center velocity, Eq. (|25|) 
below, one needs the following quantities: 

R(0) r (0) r ( ) 

b (o) _ ^_ _ £^ efl , f£l e _ 6 (o) , 6 (o) (6) 
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Af) = A(°) + ^||b( ), (7) 

,m = , (R (o) , (?aA„2 



e v <j>T = ^ + W jf, (8) 
E (o) x B (o) 

= c 3— = 0, 9 

Ef ) = -VC (0) = -- (5 (0) )'e, (10) 



and 



with 



(h {0) Y 

Bf > = V x AJ<« = rfb + ^.ill- (e r x b<°>) , (11) 

^11 p V ' 



and 



£Q 0) = B*<°> ■ b<°> = + ^ W |,y^ (12) 

Y ez (r) = b<°> • (V x b<°>) = (6S 0) )'6f - (b^)'b^ + (13) 
With the aid of Eqs. (|6|-p~3|) the guiding center velocity takes the form 

- ^ m " T^W^ (* >< b< °>) + ^ (* >< b,0> ) ■ (") 



to,. r 



e u Blf 

with cu* = cr ^" . The first, second and third terms in ([14]) are the component 
of v^J parallel to B^°\ the grad-B drift and the curvature drift, v^j) has no r- 
component and therefore r is a constant of motion. Since there is also no force 
parallel to B^, another constant of motion is the parallel guiding center velocity 
v\\. The guiding center distribution functions f^J are therefore functions of r, v\\ 
and the magnetic moment /i. 

To calculate the current density J^ -*, we apply the general formula (8.15) of 
Ref. ||18|| ) which was derived in the context of collisionless Maxwell-drift kinetic 
theory. The result is 

j(o) = i V xB(°) 

47T 

V J 

- £ cV x / dv^ [BlfffJ („b - ^||V g „ ± ) } , (15) 
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where v gv ± = v gu —v\\b. The first and second sums in ( |T5D represent, respectively, 
the guiding center and the magnetization contributions to J^ -*. Taking the cross 
product of Eq. fllBD with B^ -* and using Ampere's low on the left-hand side of the 
resulting equation we obtain after some straightforward algebraic manipulations 



d_ 
dr 



p(0) + 



B(o) 



+ 



(B, 



(0)x 2 



47rr 



+ n(r) = 0, 



with 



and 



P<°> = £ / dv^ ^B^ffJ 



(16) 



(17) 



n(r) 



v r 
-2E 



(o) 
y 



m,,c 



—m u v\ 



Relation ( Jl~6| ) can also be derived by the momentum-conservation relation Tff = 
with the tensor given in explicit form by Eq. (76) of Ref. [0. (The comma 
in the subscript denotes covariant derivative.) For Maxwellian distribution func- 
tions it holds that TI^ = 0, and Eq. (|T6| ) reduces to the known MHD equilibrium 
relation. 



III. SECOND-ORDER PERTURBATION ENERGY 



The second-order energy of perturbations around an equilibrium state is given by 



d 3 x T, 



(2)0 



i 



(19) 



where Tq^° is the energy component of the second-order energy-momentum 
tensor [El 



T (2)A 
P 



' (o) (dSjp e d*H® 
U { dq« c K I dP x dP K 



h/ (Q) F (i) f^ 0) | ?W 



dP x dF T(7 



dPi 



4 



-2^ } E / dUP 

v J 



f(0) (dS?) e IPHW 
Iv \dq« c K )dP K dFl$ 

4tt 



+5 P A (E / d4dPf2>\H? - H^) + ^FWF^j . (20) 

Here, the superscripts (0), (1) and (2), respectively, denote equilibrium first- 
and second-order quantities; A p = (—0, A), where is the scalar potential and 
A the vector potential of the electromagnetic field; F^ is the electromagnetic 
tensor; are generating functions associated with the perturbations; the scalar 

quantity \f^~Qp~) results from the contraction in the second-order tensor 
Iff?' j ; the rest of the notation is defined on page 273 of Ref. [fj. 
In expression (|20]) the time derivatives ^ are given by 

_ e A m - _ ls m H m + ^ a( d . ^ _ F m??t (2 i) 

dt v ~ c <9P mA aF (0) ' 

where the mixed variable Poisson bracket is defined as 

da db da db 



[a A 



dqi dp dPi dqi 



The Hamiltonian for the guiding center motion of particle species v is obtained 
from the Lagrangian 

L v = (^) At • x - e v <j>l (22) 

with 

A* = A + g 4 b, 

e„ 



2 



e,0t = + ^5 + ((g 4 ) 2 +v Ej . 



2 

E x B 



w E = c- 



B 



2 



E = -V0--— , B = VxA, b = -. 

c at B 

This Lagrangian is defined in terms of the variables 

t, x = x (g 1 , q 2 , q 3 ^j and q 4 . 

5 



Here, q l , q 2 , q 3 are generalized coordinates in normal space and q 4 is an additional 
independent variable for wich one of the Lagrangian equations yields the relation 
q 4 = v ■ b = v\\. The momenta canonically conjugated to x and q 4 follow from 
(0) as 

p = ~d± = lkf e = T Pt = W = °' <23) 
where e' are the reciprocal base vectors. Since Eqs. (|23|) do not contain x and 
q 4 , they are constraints between the momenta and the coordinates. It therefore 
follows that Hamilton's equations based on the usual Hamiltonian corresponding 
to the above non-standard Lagrangian are not the equations of motion. To over- 
come this difficulty, Dirac's theory of constrained dynamics |^0] is applied, which 
yields the Dirac Hamiltonians: 

H v = e v <j>* v + Vg,, ■ (p - (e v /c)A*) + V 4 Pi , (24) 

from which 

x = v = ^ = v 9 „ (W) = |J-B* + ^-K x b (25) 

q* = ^ = V 4 (t, x, q 4 ) = ^^rK ■ B/ (26) 
oPa v 1 vrt v B\ 

follow. (Here, E* = V<# - B„* = Vx A* and B* H = B* • b.) Special 

solutions of the equations of motion following from the Hamiltonians (|4]) are the 
constraints (P3|). The distribution functions /^(x, g 4 , p,p4, t) must guarantee that 
these constraints are satisfied. As concerns this requirement, it is important to 
note that p — (c) A* = and p 4 = do not represent special values of some 
constants of motion. Therefore, ^-functions of the constraints are not constants 
of motion either. On the other hand, f v must be proportional to such ^-functions 
and, at the same time, also a constant of motion. Both conditions are uniquely 
satisfied by 

/„ = S( P4 )5 (p - ^A*) Bt [{ f gu (x, g 4 ,//,t) , (27) 

where the guiding center distribution functions f gu are constants of motion and 
solutions of the drift kinetic differential equations 

B -w + ^ +vt w =0 - (28) 



and 
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In the present paper, the second-order perturbation energy is calculated for 
the case of the equilibria defined in Sec. II and for initial perturbations = 
A^ 1 ) = 0. It is also shown a posteriori that one can choose initial perturbations 
without changing the particle contribution to the energy in a way such that the 
corresponding charge density p^ 1 ' vanishes. Therefore, choosing initial perturba- 
tions of this kind, we can put from the outset 



Equation (|21|) then reduces to 



dt 



0, A<P = 0. 



(29) 



and the Dirac Hamiltonians to 



#(o) 



+ V 



(0) 



(30) 



(31) 



The second-order perturbation energy (Eq. (|19D) takes then the form 

+ E f d 3 xdq*dPfM (n^-n^] 



dPi J 



with 



dspr 

OP 



_d 
dq 



f \ Ju dP J q lJu dP x 



(32) 



(33) 



x(g 1 ,g 2 ,g 3 ) = x(r, 6, z) and d 3 x = q 1 dq l dq 2 dq 3 = rdrdOdz. After a lengthy 
derivation, which is presented in Appendix A, Eq. (|3~2"D can be cast in the concise 
form 



i?(2) 



E / S{r)drdv\\dp 



_B*(o) 



/ 




,(o) 

9^ 



1 0/g> 



dvw ^ iot^ dr 



(34) 



Here, S(r) is a normalization surface (Eq. ( |A.2S| )), G^'(r, g 4 ,/i) are arbitrary 
first-order functions related to the perturbations (Eq. ( |A.25| )); k# 2 , k\\ and k± 
are the wave vector lying in magnetic surfaces (Eq. (|A.26 )) and its components 
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parallel and perpendicular to B^. We note that depends on only via 
|G«| 2 . 

Since the first-order charge density p^ is a vn and ji integral over an expression 
that is linear in and therefore also linear in G\p, one can satisfy the relation 
p^> = by a proper distribution of positive and negative values of on which 
does not depend. 

For a vanishing field line curvature (B^ = or r — ► oo), Eq. (^) reduces 
to the F^ 2) expression for plane equilibria which was derived previously 0] (Eq. 
(82) therein). New terms here are the curvature-drift component of v^y, and 



M 0) ) <9/ (0) . . <9/ (0) 



k± J( Q j r — The latter term signifies that ^ plays a role for pertur- 

(d v II II 

bations propagating not parallel to H(°\h± ^ 0), a property arising from the fact 
that the curvature drift component of v^J depends (quadratically) on the parallel 
velocity v\\. 

IV. CONDITIONS FOR THE EXISTENCE 
OF NEGATIVE-ENERGY PERTURBATIONS 

First it is again noted that the conditions for the existence of negative-energy 
perturbations hold if the chosen frame of reference is that of minimum energy. 
Perturbations propagating parallel, obliquely and perpendicularly to B 1 - ^ are 
separately considered. 

A. Parallel modes (k± = 0) 

In this case Eq. ( pip reduces to 

F (2) = -SJ2J rdrdv\\dp 



D*(0) 



GV 



(35) 



dv 

Thus, one obtains F^ < if 

Jgu >0 (36) 



dv\\ 

holds for some r, v\\ and p for any particle species v. Condition (0), first derived 
by Pfirsch and Morrison || for a homogeneous, magnetized plasma, guarantees 
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the existence of negative-energy perturbations without any restrictions on the 
magnitude or orientation of the wave vector other than k» 7^ 0: it suffices to 

localize G\P to the region in r, v\\ and /i where v\\ ^ > 0. Outside this region 

vanishes. All the other i.e. with A 7^ u, are set equal to zero. The 
sign of is then determined only by the sign of the integrand in the region 
of localization. This result agrees with those obtained by Correa-Restrepo and 
Pfirsch for several Vlasov-Maxwell equilibria | 



B. Oblique modes (kn 7^ and kj_ 7^ 0) 



With the definitions 



h 11c dB® V]] (bf)'' 
C = v\\t jl ^—77^— h 



and 



D 



h 91 



k ± dv l{ ^(°) 
Eq. (H) yields < if 



V\\ 



Or 



r dv\\ 



C > and D > 



or 



C < and D < 0. 
The following two cases are now considered: 



(37) 

(38) 

(39) 
(40) 



a) Let us first assume that 



OV\\ 



(41) 



again holds locally in r, vn and \i for any particle species v. It then follows from 
inequalities fl3~9|) and ( [40"D that 



fcii k\\ 

-!L < min(A I/ , Af„) or -!L > max(A^, M v ), 



(42) 



with 



A, 



1 /ic d5( ) 1 (4 0) ) ; 



v ll e. 
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and 

M _ vj {"f } f i (ejgy 1 

" ~ u,t m r + Jfl dr \ dg* ) ' 
The perturbations G u 1 ^ are localized as in the previous case of parallel propaga- 
tion. The orders of magnitude of A u and M v depend on the particle energy. For 
thermal particles, these being the most representative particles, it holds that 

\ kv \ « \M V \ « « 1 

(r ) 

I^Kj^jth. j s thermal Larmor radius), and consequently condition (|42f) imposes 
no essential restriction on the magnitude or orientation of kg z associated with 
negative-energy perturbations. 



b) On the other hand, if 

^<0, (43) 



dv\\ 

holds at some r, v» and [i for any u, a condition which is more frequently satisfied 
(e.g. in the case of Maxwellian distribution functions), it follows from inequalities 
and Q4TD that negative-energy perturbations exist if, in addition to 



fcii 

minfAj,, M v ) < -± < max(A v , M„) (44) 

holds. For thermal particles the latter condition implies that 

i « « L (45 ) 

k± r 

Therefore, the most important negative-energy perturbations, in the sense that 
the less restrictive condition (|43| ) is involved, concern nearly perpendicular modes. 



C. Perpendicular modes (k\\ = 0) 

In this case, with the aid of the equilibrium condition (pf), Eq. (|34]) reduces to 

TXT / 7 \ 2 

F<*> = ^S/^^l^f-^^ (J^) (46) 

with 



10 



and 

Q = (dfl _ (b¥>)>dfv\ ^ 

v \ Or r dv\\ J 

Here, W v \\ and W u ± are the parallel and perpendicular particle energies. Negative- 
energy perturbations exist whenever either of the conditions 

Ru<0 and Q v > (49) 



or 

R v > and Q v < (50) 
hold. Condition (|50"D , which cannot be satisfied by plane equilibria with singly 

jp(0) 

peaked pressure profiles for which = ^ < 0, determines a new regime of 
negative-energy perturbations. The consequences of (|49D and ( |50|) for straight 
tokamak and reversed-field pinch equilibria are examined in Sec. V. To simplify 
the notation, the superscript (0) will be suppressed on the understanding that all 
quantities pertain to equilibrium. 



V. PERPENDICULAR NEGATIVE-ENERGY PERTURBATIONS 
IN EQUILIBRIA OF MAGNETIC CONFINEMENT SYSTEMS 



A. Straight tokamak equilibria 



Straight tokamak plasmas which are close to thermal equilibrium can be de- 
scribed by shifted Maxwellian distribution functions 



f9V ~ ' 27T 



m. 



2 N u (r) 
Tu ,2 (r) 



cxp 



fj,B(r) 



l 

2m„ 



'II 



VJr) 



T v (r) 



(51) 



where N v and T u are, respectively, the number density and temperature (in energy 
units) for particles of species v. The shift velocity V p satisfies 



V v 



( r Lu, 



tli 



< 1 



and, as shown later, leads to a net "toroidal" current. 



(52) 



In the remainder of the paper the analysis will be carried out up to zeroth 

' (r Ll/ ) th 



order in (r^th/ro, i-e. small terms of the order of 



(with n > 1) will 
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be dropped. In this context, from Eq. (|1"8| ) one obtains EL, w 0, and Eqs. (|T6| ) 
and ([47]) reduce, respectively, to 

B 2 \ Bl 



(53) 

Ul \ Oil J t'/l I 

For distribution functions ([51]), negative-energy perturbations exist if the relation 

RuQu = Ru (^f\ U v f gv < (55) 



4irW v ±N v 
+ W^N> (56) 



SlniV (57) 
It is now assumed that both the density and temperature profiles are singly 
peaked and therefore r\ v > for all v. Negative-energy perturbations thus exist 
in the following two regimes: 

a) R u < . This implies that R v \ ~sf- \ > and, consequently, condition ( |55| ) is 




satisfied if U v < 0. Since the last two terms of U v are non-negative and vanish 
for W v \\ = W u ± = 0, the condition U v < can be satisfied if 

V»>1 (58) 

holds for some particle species v. The existence of perpendicular negative-energy 
perturbations for any perpendicular wave number is therefore related to the 



threshold value of 2/3 of the quantity r\ v . As discussed in Ref. |L4|], this threshold 
value is subcritical in the sense that it is lower than the critical value r\1 ~ 1 for 
linear stability of temperature-gradient-driven modes. 



b) R v > . Condition flog) is now satisfied if U v > 0. In this case nega- 



tive energy perturbations exist for any k± without restriction on the values of rj u . 
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We now find the part of the velocity space occupied by particles associated 
with negative-energy perturbations (active particles). The particular particles 
with energy components W v \\ = T u /2 and W u ± = T u , and consequently with 

velocities equal to the root mean square velocity (v u ) rms = y are first exam- 
ined. For these particles, henceforth called representative particles, the quantity 
U v becomes independent ofrj u . Condition R u < 0, U v < is then impossible and 
condition R v > 0, U v > , concerning the new regime, takes the simpler form 

4?r N v ( , Bl \ 
" 1< ^ivf( P+ 2^) <0 - < 59) 

Condition d59|) guarantees that the representative particles are active particles. 
For particles with arbitrary velocities the part of the velocity space occupied by 
active particles is determined on the basis of analytic solutions constructed in the 
following way: 



Inserting the distribution function (51) into the equilibrium equation (15) and 



carrying out the integrations with respect to vn and //, one obtains 

CO c 

Je = b e J2 e uN u V u + -^P' = -—B' z (60) 

v B 47T 

and 

- J z = -h £ e v N v V v + °^P' = -^r~ (rB e )' . (61) 

To get some simple kind of insight, we now restrict discussion to T« = 0, a case 
often considered in the literature, e.g. |nj, [^2j. For cold ions Eqs. ( |60| ) and (|61"| ) 
yield 

c^P' - eb e N e V e = ~& z (62) 

D 47T 

c b Ap> + ebzNeVe = -±\rB e )' (63) 
B Air r 



and 



with e P = —e and 



N e T e . (64) 



Let us briefly discuss here the meaning of V e : For V e = and a constant 

"toroidal" magnetic field B z = Bo one obtains from Eq. ([H]) the "toroidal" 
current density 

J z = ~%P'- (65) 
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On the other hand, Eq. ( |62] ) for this case yields P' — 0. Hence, there is neither a 
pressure gradient nor a toroidal current. For an r-dependent toroidal magnetic- 
field component, B z {r), and V e = 0, Eqs. ( |60|) and (|6l"l) reduce to 



- s (rB,) ' = l p ' (67) 

For 7^ 0, one can readily show that their solutions satisfy the relation Bg = — ^ 
(with c= const.) and therefore they are singular at r = 0. For Bg = 0, Eq. ( |67D is 
trivially satisfied and Eq. (|66|) describes a shearless stellaratorlike configuration 



with vanishing toroidal current, a case which was studied in Ref. [M 



To obtain analytic straight tokamak equilibria, it is convenient to use, instead 
of Eqs. (H) and Eq. (0) and 



which is equivalent to the equilibrium condition (|53|) . Here, ip(r) is the usual 
poloidal flux function. Assigning the V'-dependence of the functionals P(ip) and 
B z (ijj) and the r-dependence of V e (r), one obtains from the solution of Eq. (|68| ) 
the poloidal magnetic field B# = Vip xe z = — ^" e e, the electron density from 
Eq. ( |62|) and the electron temperature from Eq. (|58|). We have considered two 



classes of equilibria: 

i) B\ and P are linear in ip and ii) S 2 = constant and P = quadratic in 
ip. For both classes we chose r] e = 0, r\ e = 1, r\ e — > oo, and ?7 e < 0, the latter 
with singly peaked density and hollow temperature profiles or with singly peaked 
temperature and hollow density profiles. From these equilibria the following 
results are deduced (Two examples are discussed in Appendix B.): 



1. A substantial fraction of the thermal electrons are active, e.g.: 

• For linearly (marginally) stable equlilibria of a strongly diamagnetic 
plasma with T] e = 1, more than one-third of the thermal electrons are 
active. 

• For linearly stable equilibria of a paramagnetic plasma with flat elec- 
tron temperature profiles, the entire velocity space is occupied by ac- 
tive electrons. 
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2. The fraction of active particles increases from the center to the plasma edge. 

3. The fraction of active particles in a paramagnetic plasma is higher than in 
a diamagnetic plasma with the same pressure profile. 



B. Reversed-field pinch equilibria 



The same distribution function ([51]) is employed to derive force-free equilibria. 
Linearizing Eq. (|68| ) by means of the ansatz P' — and B z oc tp and then solving 
the resulting equation, one obtains B z = B z (0)J (p) and Bg = B z (0)Ji(p). These 
profiles satisfactorily describe the central region of the relaxed state of a reversed- 
field pinch [p3| . We note that perpendicular negative-energy perturbations do not 



exist in force-free plane equilibria with sheared magnetic field, which were studied 



in Ref . |14[] , because for this case the second-order perturbation energy vanishes. 
For cold ions and by appropriately assigning the mean electron velocity profile, 
one can derive equilibria with various density and temperature profiles having 
non-positive values of r] e for which negative-energy perturbations exist and a 
substantial fraction of active, thermal electrons are involved. 

As an example we considered an equilibrium with peaked density and hollow 
temperature electron profiles: 

K = const., N e = N e (0)^-r, T e = T e (0)^. (69) 



Condition (^) then yields 



W e± W e || 5 



for any p, which implies that more than half of the thermal electrons throughout 
the poloidal cross-section are active. 



VI. CONCLUSIONS 



The general expression for the second-order perturbation energy, derived by 
Pfirsch and Morrison in the framework of linearized collisionless Maxwell- drift 
kinetic theory, was evaluated for the case of circularly cylindrical equilibria and 
vanishing initial field perturbations. From this expression we obtained the fol- 
lowing conditions for the existence of negative-energy perturbations, which need 
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only be satisfied locally in r, v\\ and ji and are valid in the reference frame of 
minimum equilibrium energy: 



1. If the equilibrium guiding center distribution function fffl of any species v 

has the property v\\ v 9 " > 0, parallel and oblique negative-energy pertur- 
bations (fcii 7^ 0) exist with non essential restriction on k. 

9/(0) 

2. If t>|| v 9 " < 0, the oblique negative-energy perturbations possible are nearly 

perpendicular. With the quantities R v and Q v defined by (fT7|) and PS|) , the 

condition for perpendicular perturbations is R V Q V < 0. From this condition 

it follows that the curvature, which is associated with B-f \ modifies the 

jp(o) 9/(°) 

plane-equilibrium condition a \ T < 0. 

ar or 



For the case of tokamak equilibria there are two regimes: 



1. If < 0, the existence of negative-energy perturbations is related to the 

d In T 

threshold value of 2/3 of the quantity r\ v = ^\ n j\f ■ 

2. If R v > 0, a new regime appears, not present in plane equilibria, in which 
negative-energy perturbations exist for any value of r\ v . 



For various tokamak cold-ion equilibria with negative and non-negative values of 
r] e , a substantial fraction of the thermal electrons are associated with negative- 
energy perturbations (active particles). In particular: 



1. For linearly (marginally) stable equilibria of a strongly diamagnetic plasma 
with i] e = 1, more than one-third of the thermal electrons are active. 

2. For linearly stable equilibria of a paramagnetic plasma with flat electron 
temperature profiles, the entire velocity space is occupied by active elec- 
trons. 



The part of velocity space occupied by active particles increases from the center 
to the plasma edge region and is larger in a paramagnetic plasma than in a 
diamagnetic plasma with the same density and temperature profiles. 
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It is also shown that, unlike in plane equilibria, negative-energy perturbations 
exist in force-free, reversed-field pinch equilibria with a substantial fraction of 
active particles. The present results, in particular the fact that a threshold value 
of r\ v is not necessary for the existence of negative-energy perturbations, enhance 
even more the relevance of these modes. 
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APPENDIX A: DERIVATION 
OF THE SECOND-ORDER PERTURBATION ENERGY 
FOR CIRCULARLY CYLINDRICAL EQUILIBRIA (Eq. (|3j)) 



We start from the expression d32|) . In order that the constraints ( p3|) be satisfied, 
the terms in the first sum of Eq. (|32|) , which contain derivatives of fj)°\ are 
integrated by parts: 



/ d 3 xdq 4 dP— v —— = d 3 xdq 4 dPfP — ^t^^t 

J H dt dq l V dPi ) J dPi dq l dt 

+ J dqWdq'dPf^^t^ 



Furthermore, because of 



Eq. (10) of Ref. yields 



d 2 nw 

3P l dP k 



(A.l) 

0, (A.2) 



H {2) = 0. (A.3) 

Using Eqs. flSOJ), (f53|), ( jOQ , (|A~^ ) and Eq. (12) of Ref. (H| for H^ 2 \ and 
noting that the contribution of the last term in ( |33"D cancels the contribution to 
of the last term in ( |A.1| ), Eq. (|32|) is put in the form 

F{2) = E / d?xdq 4 dPfl 0) A, (A.4) 
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with 



A 



_ dS^d_ ( dH® dSW \ _ 1 
dPi d$ V <9<f dPj J 2 
dffl d ( dSW dH™ \ 
dPi dq i V dqi dPj J 



asp asp 

' dPi dPj 



(A.5) 



( i,j = 1, 
vector 



4). To make treatment of the constraints easier, we introduce the 



V = — 



m,. 



^ Af ) (x, g 4 ) 



It can then be shown that 



dqi 



0. 



v=o 



H (o) 



■ 1.1 



d 2 H^ 



v=o 



v=o 



and, consequently, 

(^ 0) ) 



dqi dqi 

dSjpdSW dSW d fdHjpdSp 



w dPi dP j dPi dqi \ dqi dP j 



(A.6) 

(A.7) 
(A.8) 

(A.9) 



We note that the constraint P4 = is not involved here, because P4 does not 
appear in (Eq. (0)). With the aid of Eq. (pOl ), Eq. (TAJ) is now written 

as 



ldS® d (dH^dS^\ dSjp d (dSjpdHj®\ 



dPi dqi 



{ dqi dPj J 



2 dPi dqi \ dqi dPj 
The two terms of Eq. ( |A.10| ) will be calculated separately. 

The first term can be written as 

dS^ d ( dHjpdSpA _ d 2 H^ dS?> dSW 
dPi dqi \ dqi dPj ) ~ dqidqi dP { dPj 
d 2 H^ dS<» dSP 



(A.10) 



d (dH® \ dS<pdS<p 



QqkQqi 



v v _ . v 1 v v 

dP k dPi W I dP l ) dq l dP< 



d(q 4 ) 2 { dP 4 J 



(A.ll) 



4 and k, I — 1, . . . 3. Since the equilibrium quantities depend just 

d 2 H<® , d (dH<®\ 



with i,j = 1, 

on g 1 , the only non- vanishing components of 
to Hamiltonians (|3"TD, are 

a 2 H® 



WW and w{~^t)' accoT(hns 



d{ q i y 



v=o 



c dq 1 dq 1 



(A.12) 



18 



and 



v=o 



_d_ ( dH® \ 

dq 4 \ dq 1 J 
On the basis of relations (|A.12|) , (|A.13|) and 

d 2 H® 



d{A 



dq 1 dq 4 



d{A 



v=o 



(dq 4 ) 2 
Eq. (| A. 11| ) reduces to 

dPi dq 1 \ dqi dPj J 



c dq 4 dq 4 



(A.13) 



(A.14) 



d{v$)) 1 d(A$°)] 
dq 1 dq 1 



d(vmd(A<% e v d{vf)) l d{Af» 



c dq 1 dq 4 



c dq 4 dq 4 



4 



We now calculate the second term of Eq. (|A.10|) . By virtue of —gft 
second term on the right-hand side of 

dS^ d f dSP dH® \ 
dPi dq* \ dqi dPj J 



(A.15) 



0, the 



dPi dq i \ dq 1 dPi ) 
, dSV d (dSV>dH®\ 



dP, dq* V dq 4 <9P 4 ' ^' 16 ^ 
= 1, . . .4, I = 1, ... 3) vanishes. We note here that, whereas Eq. ( |27| ) for 
f v is sufficient in the nonlinear theory to pick out the correct solutions, this is 
not so with the linearized theory. In this case, since the constraints are imposed 
along the perturbed orbits, a displacement vector (£,£4) in x, q A space, similar 
to that in macroscopic theory, is introduced 0; that is, since the zeroth-order 
distribution function always selects V = and P4 = 0, with V as defined by Eq. 
( |A.6| ), it is reasonable to expand in powers of V and P4. 



S« (x, g 4 )-$.m„V-^P4 
+ higher-order terms, 



so that 



dP 



V=0, P 4 =0 

Using equation ( [A. 171) , one has 



dq 1 



asp 



dq 1 

asp 

dq 1 



dP d 



dP k 



(A.17) 



(A.18) 



dq 



V=0, P 4 =0 



asp 



OP, 



e v d(A*W) k 



c 



dq' 



asp 

dP k 



(A.19) 
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and, therefore, 



dq 1 



OP, 



dq 1 



,/o)y 



d(A*W) k DSP 



v 



dP k 



. (A.20) 



Since A*(°) depends only on q 1 and v^j is perpendicular to e r , the last term 
on the right-hand side of Eq. ( |A.20|) vanishes. This has the consequence that 
higher-order terms in expansion ( |A.17| ), after the constraint V = is imposed, 
do not contribute to Eq. ( A.22|) below. Applying the operator (m = 1,4) 

to Eq. QOOD , one has 



d 



( dSP 



p <9P, J 



dq m \ dq l 

e v d{A<°)) k d 



dq r ' 



_d_ 

dq n 



_d_ 

dq n 



(o)y 

qv I 



dP k 



dq l 



dq r ' 



(o)y 

QV I 



as?" 





d 




dq m 




(o)y 

gu ) 


V- 


V 


'* ( 


dP k \ 



(o)y 



IV 



vJ 



VJ 



V/ 



dq m 



t as® 




\ dq k 


J 



With the aid of expansion ( |A.17| ) and (|A.18 ), the last equation yields 



_d_ 

dq n 



dq 1 



<9P, 



v=o 



dq" 



9^ 



c <9g m 

Equation (|A.16|) can then be written in the form 

dPi \ dqi dPj ) " dq 1 dq 1 
K 9V> dq k dq iq dq 4 dq 1 S 



<9g 



(A.21) 



(A.22) 
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K 9U> d<fdq 1 ^ c K 9u) dq 1 <9g< 

+ c (V) dq* dq^' [AM) 
Inserting of Eqs. flQ5D and (TQ5|) into Eq. (|O0|) leads to 



2 c dg 1 dq 1 ^ ' c dq 1 dq 



le 



2 c <9g 4 <9g 4 ^ ' dq 1 dq 

<9g fc <V <9g 4 <9g' ? 

A^ffi^i + f* < 4 (A 24) 

with k, I = 1, . . . 3. 

Since the equilibrium is independent on g 2 and q 3 , an appropriate ansatz for 
the functions is 

Si 1] = G[ / 1) (g 1 ,g 4 ,/i)e i(k23x) . (A.25) 

The wave vector k 23 = kg z introduced here has constant covariant components 
&2 and k 3 and physical components kg and 

k 2 3 = ^2 7^ + k 3g^ = k S e + k z^z = Kz- (A. 26) 

Therefore, it lies in magnetic surfaces. We rewrite the integral over the momen- 



tum space according to the rule [18 



and introduce real quantities by 

AB -> ~3M*fl; (A.27) 

then inserting Eqs. (37-40) of Ref. for «f and Eq. (|OE| ) into Eq. ( [04]) , 
integrating with respect to g 2 between g 2 . and g^ + ^r- an d with respect to g 3 
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between q$ and q$ + -77-, taking into account that d 3 x = q 1 dq 1 dq 2 dq 3 and denning 
the normalization surface S^g 1 ) by the relation 



Siq 1 ) = q 



f q%+2n/k 2 f4+2n/k 3 

1 1 dq dq 6 , 



(A.28) 



Eq. (|A.4j ) (after a lengthy algebra) can be written as 



c 
2^ 



9 



with 



5 



23 



dq 



B 



(A.29) 



9^ 
I dg 1 



1 c /avg gAg) ' 

2 e|/ B*,| 1 dg 1 ' dq 1 



1 -L 



fc 2 - 



(if ' k -> 



2m yJ B 



HI 



2:-! 



^ . b <»> 

,4 



^ dq 4 



^( B » <0, -M(^-^)}| G «| 2 /<?, 

M*? 1 ) = (b<°> x k 23 ) • = (b< > x k, 2 ) • e r 



(A.30) 
(A.31) 



and 



A;,, (g 1 ) = b<°> ■ k 23 = b 2 k 2 + b 3 k 3 = b e k e + b z k z . (A.32) 

Integration by parts of the terms in (|A.29 ) in which d\G^\ 2 / dq 1 and d\G^\ 2 /dq 1 
appear, leads to the expression 

2 ®fgu 



F (2) = - E / Stf)<Pxdfd»yL (vj ■ k 23 ) (B^ ■ k 23 ) |G« 



2c 



(vg? • k 23 ) fc ± |c#> 



2 ^ + # + C 



dg 1 



dq 4 
(A.33) 



with 



/9v(o) 



2e u [ Oq 1 ' k 



23 I 1 H — — fv 
7 2e„V 



2:-! 



dg 1 



,(0) 



(A.34) 
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Inserting Eqs. (|7|), (|TT| ) and ( |i4|) for, respectively, A*(°\ B*(°) and viy into Eqs. 



( |A.30| ) and ( |A.34j ) and using the identities 

*. s _ yfcfcl + M _ 2 ^ s _ M,„ + ^ (A , 5) 

dr r r r r 



and 

^ + ^_i^ ± ,y fc% (A .36) 
one can show (after tedious but straightforwrd algebraic manipulations) that 

B + C = 0. (A.37) 
With the aid of Eq. QQ7D , Eq. ( ggg ) reduces to Eq. 



APPENDIX B: ACTIVE PARTICLES 

The part of the velocity space occupied by active particles is determined by 
means of analytic solutions of Eq. (]68f) . With the ansatz ^ (P + B 2 ) = const., 
the solution of Eq. (|68|) is of the form ip oc p 2 , with p = r/r . This yields a class 
of equilibria with the following characteristics: 
Peaked parabolic pressure profile 

P = P(0)(l-p 2 ), 

where a is a parameter which describes the magnetic properties of the plasma, 
i.e. the plasma is diamagnetic for a 2 < 1 and paramagnetic for a 2 > 1; 



B z 



P z 2 (0) + 8vrP(0)(l-aVl 1/2 ; (B.l) 



B e = 2^7rP(0)ap; (B.2) 

constant "toroidal" current density. Assigning appropriately the shift electron 
velocity profile, one can construct equilibria with a variety of values of i] e . Two 
examples are discussed below. 

1. rj e = 1 equilibrium 

Choosing the V e profile as 

Ve = v e (0)^ f (l-P 2 r 1/2 , (B.3) 
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with 



and 



B 2 = (B 2 + B' 2 Z ) = B 2 (0) + 4ttP(0)(2 - a 2 )p 2 

1/2 



B 2 Z {0) + 4vrP(0)(2-a 2 )p 2 



one obtains 



N e = N e (0)(l - p 2 ) 1 / 2 and T e = T e (0)(l-p 2 ) 1 



/2 



(B.4) 
(B.5) 

(B.6) 

-1/2 



Therefore, r] e = 1 holds for all p. We note here that, owing to the (1 — p 2 
dependence of V e , the equilibrium profiles are possible only in the interval < 
P < Ps < 1) with p s appropriately chosen so that inequality fl52|) is satisfied (e.g. 



Ps 



with 



. Condition (|59|) , concerning the representative particles, yields 
< (3(1 - p 2 )(a 2 - 1) < 1, 



P(0) 



P(0) 



const. 



(B.7) 



(B.8) 



]^)) must 
Thus, the 



P 2 /8tt P 2 (0)/8tt 
The requirement that the "toroidal" magnetic field modulus (Eq. 
be non negative sets the upper limit 1 + (3~ x on the values of a 2 
right hand inequality of condition ( |B.7| ) [(3{1 — p 2 )(a 2 — 1) < 1] is satisfied for all 
possible values of a 2 . The left hand inequality [0 < (3(1 — p 2 )(a 2 — 1)] is satisfied 
for a 2 > 1 and, therefore, only in a paramagnetic plasma the representative 
particles are active. For particles with arbitrary velocities, conditions ( ff5| ) and 
(pOD, respectively, yield 



W, 



e± 



< - ( — - 1 ) and 



1 - -(3(1 - p 2 )(a 2 - 2) 



and 

W 4 1 
W^ > 2 



a' 



1 ) and 



1 



(3(l-p 2 )(a 2 -2) 



P 



W P 



l-(3a 2 (l - p 2 



P 



1 - (3a 2 (l - p 1 



1 

Pe 2 

B.9) 

W«n 1 
> -. 



P 
(B.10) 

For a strongly diamagnetic plasma (a — > 0) condition ( p.lOj) is impossible and 
condition (|B.9|) yields 



l + (3(l-p 2 



— — H ^ < -. 

T, T, 2 



(B.11J 



The part of the velocity space occupied by active particles is depicted in Fig. 
1. We note here that in this and in the following figures the dotted area stands 
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for the active particles at plasma center (p = 0), while the area filled by circles 
stands for the additional part of active particles at p = p s . 




Figure 1: The part of the velocity space occupied by active electrons for a 
strongly diamagnetic plasma with r\ e = 1, wich is deduced from Eq. ( |B.11| ) 

Mp) = i + /3(i- p 2 )]. 



Since (3 is one order of magnitude lower than unity, relation ( B.11J) implies that 
nearly one third of the thermal electrons are active. In addition, the fraction of 
active electrons slightly increases as one proceeds from the center to the edge, 
because the factor 1 + (5(1 — p 2 ) which multiplies W e ±/T e in relation ( p.ll|) is a 
decreasing function of p. 



For an equilibrium with constant "toroidal" magnetic field (a 2 = 1) conditions 
flBDD and flBlCl) reduce to 



< - and 

W e± 2 



and 



1 ! ' > - and 



W e± 



1 + \P(1-P 2 ) 



W, 



We_ 



+ 



1-/?(1-P 2 



1 

K 2 



1 

2 



(B.12) 



(B.13) 
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The fraction of active electrons, following from conditions ( |B.12| ) and ( |B.13| ), is 
depicted in Fig. 2. 




Figure 2: The part of the velocity space occupied by active electrons for the 
equilibrium with r\ e = 1 and B z =constant, which is deduced from Eqs. ( p. 12]) 
anddFTTl [ ai (p) = 1 + 1/2(3(1- p 2 ), b x (p) = 1 - - p 2 )]. 



Nearly half of the velocity space is now occupied by active electrons. In addition, 
one can readily show that active electrons increase from the center to the edge. 

For a diamagnetic plasma with a 2 = 2 condition ( |B.9|) is impossible and 
condition (|B.10| ) leads to 



We 



+ 



1 - 2(3(1 - p 2 ) 



w n l 

T P > 2' 



(B.14) 



The fraction of active electrons is depicted in Fig. 3. 
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Figure 3: The part of the velocity space occupied by active electrons for the 
equilibrium of a paramagnetic plasma with r] e = 1, which is deduced from Eq. 

(|KTJ [b 2 (p) = 1 -2(3(1 -p% 



Nearly two thirds of the velocity space is now occupied by active electrons. As in 
the cases of a strongly diamagnetic equlibrium and an equilibrium with a constant 
'toroidal" magnetic fileld, the fraction of active particles increases from the center 
to the edge. 

2. rj e = equilibrium 



With the choice 



one obtains 



K = K(0)-||(l-p 2 )- 1 (B.15) 



N e = N e (0)(l - p 2 ) and T e = T e (0) = const. (B.16) 
Condition ( [j9|) concerning the representative particles leads to 

0<(3(a 2 -l)<l (B.17) 

and, therefore, as in the equilibrium of Appendix Bl, only in a paramagnetic 
plasma are the representative particles active. 



27 



For particles with arbitrary velocities, conditions ( |4"9"D and ( |50| ) respectively, 
yield 

— ^<- — -1 and 7 2-a 2 -r + V-^<-l B.18 
W e± 2 \a 2 J 4 V ' Te 2 ' 

and 



£l > i ( A _ A md £ (2 _ a ^ + PjVl > _ L (B.19) 
We ± 2 Va 2 / 4 V J T e 2 T e K J 

Condition (p,18| ) is, as expected, impossible for any a, because rj e takes its lowest 
non-negative value well below the subcritical one. For a — > condition ( B.19 ). 
concerning the new regime of negative-energy perturbations, is also impossible 
and therefore no negative-energy perturbations exist in a strongly diamagnetic 
plasma. For a 2 = 1 condition (|B.19|) yields 

and, therefore, half of the velocity space is occupied by active electrons for all p. 
For a 2 = 2 condition (|B.19| ) leads to 

>0 and (B.21) 



W, 



e± 



and therefore all particles are active. Thus, since the value r] e = is far lower 
than the critical value for linear stability (n c e w 1), negative-energy perturbations 
involving a large number of thermal electrons exist in a linearly stable regime. 
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